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The description of the inclusive production of hadron jets is one of the successes of 
pertubative QCD (pQCD). Quantitative agreement between data and theory has been 
achieved for jets produced over a wide kinematical range. In particular, data from CDF 
and D0 Collaborations at the Fermilab Tevatron on inclusive single jet production for 
^/s = 1800 GeV |l|, 0] are in agreement with pQCD for jet transverse energy ranging 
approximately from 15 to 400 GeV. Over this transverse energy range the cross section 
falls by seven orders of magnitude. 

With this good quantitative description of the dependence on the parameters of the 
produced jets at fixed total energy of the collision, the natural next task is to examine 
the dependence of the production cross section on the total energy. Dimensional analysis 
and scaling hypothesis predetermines this energy dependence. pQCD (see, e.g., review 
0) dictates a particular mechanism of scaling violation involving a hadronic scale, Aqcd, 
which yields a specific non-trivial energy dependence. Any deviation from this prediction 
would manifest an inadequacy of the pQCD framework for managing nonperturbative 
physics (soft hadronic radiation). 

We will show that existing data already contain evidence for additional non-pQCD 
effects which are consistent with the Balitsky-Fadin-Kuraev-Lipatov (BFKL) Pomeron 
framework. We can begin with a more general observation that a potential mechanism 
for such non-pQCD scaling violation is implied by resummation of the leading energy 
logarithms of QCD (for a recent review, see Ref. 0). In this case, the dependence of 
the cross section on the QCD running coupling constant as differs from the simple power 
dependence of pQCD. This change away from a power dependence yields an altered energy 
dependence of the cross section. 

There are preliminary data from CDF ^ on the cross section for inclusive single jet 
production at ^/s = 1800 GeV and y/s = 630 GeV. In particular, we examine the ratio 
of cross sections scaled by the jet transverse energy E_i_, taken at the same values of 
x± = 2Ej_/y/s and the rapidity of the jet rj, but at different total energies J^, |[: 
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Without scaling violation the scaled ratio is unity, regardless of the dynamics. On the 
pQCD leading order predictions for the ratio see 0. 

Comparison of the next-to-leading order (NLO) pQCD prediction |]10| with the data 



shows a noticeable discrepancy at small x± (Fig. 1). This problem has been already seen 
in previous CDF data p|. It is unlikely that complications connected with jet algorithms 



and various uncertainties of pQCD approach |TT[ may account for the discrepancy. 

In this paper we show that the discrepancy in the scaled ratio at small x± is accounted 
for by the BFKL Pomeron |Q. 

In our approach, the dependence of the inclusive jet cross section on a;^ may be ex- 
pressed as a multiple integration and summation over conformal dimensions and conformal 

-asf^{{u,n}) 



spins of the BFKL Pomerons |12|. Each term of this "sum" depends on as as Xj 



where {z/, tt,} is a set of conformal dimensions, i/j, and conformal spins, rii. To get the 
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Figure 1: The x^- dependence of the scaled cross section ratio (630 GeV)/(1800 GeV). 
Only statistical errors are shown. 



weights with which these contributions enter the "sum" as well as the n}), we use 
the effective Feynman-like rules defined in Ref. |]12[. Substituting the running 0:5 into 



the BFKL formulas for the inclusive cross sections one obtains a new x^-dependence for 
the scaled ratio R. The calculation of the scaled ratio along this line gives the result 
presented in Fig. 1. 

As was pointed out in Ref. , it is important to keep track of the most forward and 
backward jets of the events for the BFKL kinematics. In particular |T^, the inclusive 
single jet cross section is a sum of three terms. The first (second) term comes from 
the processes with untagged most forward (most backward) jet and corresponds to the 
diagram of Fig. 2a (Fig. 2b). The third term (Fig. 2c) corresponds to the processes 
with both most forward and most backward jets untagged. The analytic expressions 
corresponding to the diagrams of Fig. 2 read as follows: 
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(c) 

Figure 2: Diagrams for inclusive single jet production, x denotes tagged jet, > marks 
most forward (backward) jet. 



where are connected with the light-cone components of the produced jet momentum: 
= ±log {k^/^R) connected with the light-cone components of the momenta 

for the colliding hadrons A and B: x\ b = ± log Fa,b are the effective parton 

densities of the colliding hadrons A and B; fin is the normalization point for both the 
parton densities and the running coupling constant as (subscript H stands for "Hard"); 
fiji is the normalization point for the energy logarithms which are resummed by the BFKL 
Pomeron (subscript R stands for "Regge"). Integrations over conformal dimentions Ui ^ 
are to account for different patterns of gluon radiation occupying the rapidity intervals 
spanned by the tagged jet and the most forward (backward) jet. Note that the radiation 
involves infinite number of radiated gluons. The following analytic expressions for Wa,b, 
Wa,b may be obtained from the diagram technique of Ref. [|12| by integration over 
parameters of the untagged most forward or (and) untagged most backward jet(s): 
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^ 1 2Z/^ 

WA,B{xi,X2,i^,fiH) = W^A,i?(a;i,a;2,i^,/iH):prrr— (5) 

i (^2 + ^'^Z 



Wa,b (xi , X2 , i^, /iff) = — i 



X 
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(6) 



- Ma,b{ \-iu,HH,m{xi,X2)) 

where cc;(0, u) is the Lipatov eigenvalue uj{n, v) = {2asNc/'n')(tp{l) — Reilj{^-^^^^ + iv) 
taken at zero value of the conformal spin n {ip here is the logarithmic derivative of the 
Euler Gamma-function); 

Ma,b{\^'H,i) = I dzz^~^ FA,B{z,fj,H) (7) 

are incomplete moments of the parton densities; and m(xi,X2) = min{xi,X2}- 

One more object, Ry^{0, u), entering Eq. 4 is an element of the diagram technique of 
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+ r(M + i+^zy)' 
taken at = 0. 

The sum of Eqs. 2-4 gives the cross section {s / tt^) (f a/ dx'^dx~ which is easily con- 
nected to the invariant cross sections entering Eq. 1. 

The first calculational task is to compute the incomplete momenta of the parton 
densities entering Eq. 6. For the effective parton densities, F = g + 4:/9{q + q) [Q, we 
use the following parameterization 

Fiz,fiH) = aMz-'^'-^Kl - zr^''"\ (9) 

where we omitted the subscripts A and B, since, in our case, both colliding hadrons {p 
and p) have the same effective parton density. With this parameterization the incomplete 
moments are analytically calculable as 

M{\ij,H,l) = a{nH){B{\-h{iiH),c{iiH) + 1) 



A - fe(/i//) 



$(A - h{^iH), -cM, A - bM + 1, e-^) , (10) 



where B is the Euler Beta-function and $ is the hypergeometric function. Parameters 
a^fin), bifin), c{hh) were found by fitting the CTEQ4L ||15| parton densities. 
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The next step is to perform the integration of Eqs. 2-4 over conformal dimensions. 
This has been done numerically. The relative error of this numerical computation does 
not exceed 10% for the results presented, though, for most values of the parameters, it is 
less than few percents. 

The only free parameters which we have is the normalization points fin and fiR. They 
were taken to be proportional to the transverse energy: fiH,R = C,h,rE±, as in pQCD 
|10| . Our main result is that with fixed values of $,h = 1 and = 0.5, one has a good 
description of data p simultaneously for jet cross section at a/s = 630 and 1800 GeV 
(Fig. 3), and for the scaled ratio (630 GeV)/(1800 GeV) (Fig. 1). As is seen from Fig. 
1, pQCD fails to do this for the scaled ratio at low x±. (Note, that CDF data are still 
preliminary and require final analysis.) 
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Figure 3: Inclusive jet production cross sections for Tevatron energies. 



Qualitatively, the above result is a consequence of the energy dependence of the gluon 
radiation (see Fig. 2) which is neglected in the finite order pQCD calculations: The 
radiation is more significant for higher energies where more "diffusion" of the transverse 
momenta is allowed 

Here we should comment on the applicabilty of the running coupling constant in 
conjunction with the BFKL Pomeron and on the normalization point dependence of the 
BFKL predictions. 

There are no ultraviolet divergencies in our resummed leading energy logarithms and, 
hence, there are no running coupling constant and ultraviolet renormalization group in- 
variance. Presumably, both should arise via resummation of subleading energy logarithms 
along with appropriate new factorization theorems. Connection between "Hard" and 
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"Regge" normalization points should also be established via analysis of subleading energy 
logarithms. Our work assumes the eventual success of such a program. For recent works 



on resummation of subleading energy logarithms see Ref. [IG] and references therein. In 



addition, attempts to find the relevant versions of factorization theorems have appeared 



17| , |T8|, |T9|. For more examples of the BFKL Pomeron in phenomenological applications 
encountering the problem of normalization point dependence see Refs. [0, |19|, |2^. We 
note also that significant normalization point dependence and a deviation from the data 



2T| was observed in the calculation of forward jet production in NLO pQCD for deep 



inelastic scattering at the HERA energies, while the BFKL calculation agrees with 



the data 21 



Our conclusion is that the BFKL prediction has more capacity than the prediction of 
pQCD in fitting the data on inclusive single jet production for both ^/s = 630 GeV and 
^/s = 1800 GeV. In both approaches, BFKL and pQCD, there is a fitting parameter fi^', 
BFKL approach has additional parameter, fiji, wich was also used to fit the data. In the 
case of pQCD the good fit for ^/s = 1800 GeV is achieved at the expense of the fit for the 
lower energy at small x±. However, the BFKL prediction is able to accommodate both 
energies. 
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